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Abstract. A compartmentalized surface model of Nambu and Goto is studied on triangulated spherical 
surfaces by using the canonical Monte Carlo simulation technique. One-dimensional bending energy is 
defined on the skeletons and at the junctions, and the mechanical strength of the surface is supplied 
by the one-dimensional bending energy defined on the skeletons and junctions. The compartment size is 
characterized by the total number L' of bonds between the two-neighboring junctions and is assumed to 
have values in the range from L' = 2 to L' = 8 in the simulations, while that of the previously reported 
model is characterized by L' = 1, where all vertices of the triangulated surface are the junctions. Therefore, 
the model in this paper is considered to be an extension of the previous model in the sense that the 
previous model is obtained from the model in this paper in the limit of L' ^ 1. The model in this paper is 
identical to the Nambu-Goto surface model without curvature energies in the limit of L' oo and hence 
is expected to be ill-defined at sufficiently large L' . One remarkable result obtained in this paper is that 
the model has a well-defined smooth phase even at relatively large L' just as the previous model of L' — + 1. 
It is also remarkable that the fluctuations of surface in the smooth phase are crucially dependent on L' ; 
we can see no surface fluctuation when L' < 2, while relatively large fluctuations are seen when L' > 3. 

PACS. 64.60.-i General studies of phase transitions - 68.60.-p Physical properties of thin films, nonelec- 
tronic - 87.16.D- Membranes, bilayers, and vesicles 



1 Introduction 

One of the interesting topics in membrane physics is to 
understand why so many varieties of shapes are observed 
in biological membranes and in artificial membranes [Ij 
[SlElllj. A large number of theoretical or numerical stud- 
ies have been conducted on this problem [SMFfllSM fTUlfTTl 
[T2lll3j . Statistical mechanical view points have shed light 
on the shape transformations in membranes, and the phe- 
nomena are considered to be understood within the theory 
of phase transitions [M] . 

The well-known conventional surface model is the one 
of Helfrich and Polyakov and defined on triangulated sur- 
faces ^I ^ ITB I ITT] . The Hamiltonian is given by a linear 
combination of the Gaussian bond potential and the two- 
dimensional bending energy. The model is known to un- 
dergo a collapsing transition and a transition of surface 
fluctuations on spherical fixed-connectivity surfaces [T51 
fT9 l l2C p T 1 l22 l [23]. The smooth spherical phase and the coh 
lapsed (or crumpled) phase can be seen in the model, 
however, the above mentioned variety of phases is not ex- 
pected in the conventional surface model. 

It was recently reported that a multitude of phases can 
be seen in several surface models [24l[25l[26l[27l[28l[29l[30l . 



which are slightly different from the conventional model 
and have some additional lattice structures. The surface 
models that have a multitude of phases can be classified 
into two groups: The first is a class of Nambu-Goto sur- 
face models, where we call a surface model as Nambu-Goto 
model if the area energy term is included in the Hamil- 
tonian as the bond potential term. In fact, the Nambu- 
Goto model with intrinsic curvature energy has a variety 
of phases not only on fixed-connectivity surfaces but also 
on fluid surfaces f24l[28]. Moreover, the model with one- 
dimensional bending energy was shown to posses a rich 
variety of phases The second group of models that 
have a variety of phases is a class of compartmentalized 
fluid surfaces [301I3T]. where the surface shapes are main- 
tained by a one-dimensional bending energy defined on 
the compartments. The free diffusion of vertices due to 
the dynamical triangulations is confined inside the com- 
partments in those models. 

Therefore, the class of Nambu-Goto surface models is 
quite interesting to study their phase structure as a model 
for shape transformation phenomena in membranes. As 
mentioned above, we reported in [55] that the Nambu- 
Goto model with the one-dimensional bending energy has 
a rich variety of phases on triangulated spherical surfaces. 
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However, the model in [29j has no compartment, bet 
all of the vertices correspond to the junctions and h' 
the total number L' of bonds between the junctions is 
L' = 1. In a surface model with one-dimensional h 
ing energy, one can always assume the compartmenta 
structure. Thus, it is natural to extend the model ir □ 
by increasing the number L' from L' = 1 to non-u 
numbers L' >\ so that the model has the compartme 
ized structure. 

The notion of well-definedness for surface models s 
be noted. A triangulated surface model can be cal' 
well-defined one if any physical quantities are finite. 
Nambu-Goto model with the conventional two-dimenf 
bending energy is well-known as an ill-defined model 
cause the vertices locate in an anomalously extende 
gion in R'^ and as a consequence, the mean square 
X'^ — 'Yliii^i" becomes numerically infinite, whe 
is the center of mass of the surface. The surfaces ai 
defined also when the relation such as 5*1 = 3/2 is co 
erably violated, where 5'i = 3/2 is expected from the 
invariant property in the partition function of the su 
model. This ill-definedness comes from the fact tha 
area energy term imposes a constraint only on the 
of triangles, and hence the surfaces are dominated b 
finitely thin and infinitely long triangles [32] . 

The purpose of this paper is to see whether the 
definedness in the model in |29j is preserved or not ^ 
the compartmentalized structure is introduced. It is 
aimed at seeing the phase structure. We concentral 
the region of large bending rigidity, where the su 
shape is expected to be smooth spherical if the mo( 
well-defined. 



2 Model and Monte Carlo Technique 

On a triangulated spherical lattice, a sublattice structure 
is assumed for a discrete Hamiltonian to define the surface 
model. The triangulated lattices assumed in this paper are 
identical to those for a compartmentalized fluid surface 
model in [30] ■ We briefly figure out how to construct the 
lattice. We start with the icosahedron, which has 7V = 12 
vertices of coordination number q = b. Firstly, by dividing 
the icosahedron uniformly, we have a triangulated lattice 
of size A^= 10^^-1-2, where t is the number of partitions of 
an edge of the icosahedron. Secondly, a sublattice of size 
Nj = 10m? + 2, where m divides t, is assumed on the lat- 
tice. Then, all the edges of the sublattice are divided such 
that the sublattice has the vertices common to the origi- 
nal lattice, and this makes a compartmentalized structure. 
The total number of vertices Ns on the compartment is 
given by Ns — 30m^, and the total number of junctions 
Nj on the surface is identical to the size Nj = 10m^4-2 
of the sublattice [30]. Thus, the lattices are characterized 
by integers (^,to) or equivalently by {N, Ns, Nj) and are 
characterized also by the total number L' of bonds be- 
tween the junctions. The number L' is given by L' = L—2, 
where L = £/m— 2, which is identical to that introduced in 

m- 




Fig. 1. (Color online) Snapshots of surfaces with skeletons of 
(a) L' — 8, (b) L' — 4, and (c) L' — 2, whose corresponding 
integers {N,Ns,Nj) are (2562,882,42), (2562,1602,162), and 
(2562,2562,642). (d) A combination 1-ti • t2 is included in 
the bending energy with the weight of 1 at the q — 6 junctions, 
and (e) a combination 1 — ti ■ (t2-|-t3)/2 is included in the 
bending energy with the weight of 1/2 at the g = 5 junctions. 
The combination 1 — ti • t2 is defined also on the linear chains, 
which are connected by the junctions such as those shown in 
(d) and (e). 

Figures[Tl^a)-[IJc) show the lattices with compartments 
characterized by the spacing L' — 8, L' — 4, and L' = 
2. The corresponding integers {N, Ns, Nj) are given by 
(2562, 882, 42), (2562, 1602, 162), and (2562, 2562, 642). Thick 
lines in the figures are the compartment boundary, which 
we call skeletons. The skeletons are composed of the linear 
chains and the junctions; the linear chains are connected 
by the junctions. The vertices in the skeletons make the 
above mentioned sublattice structure. 

The Hamiltonian of the model is given by a linear 
combination of the area energy term 5*1 and the one- 
dimensional bending energy ^2 such that 

S^Si + bS2, 5i = ^A4, S2^Y.^l-U-t^), (1) 

A ij 

where b is the bending rigidity. in Si is the area of the 
triangle A, and in S2 is a unit tangential vector of the 
bond i. in 5*1 is the sum over all triangles A, and Y]jj 
in 5*2 is the sum over bonds i and j on the skeletons. We 
should emphasize that ^2 is defined only on the skeletons, 
which are the linear chains and the junctions. 

The definition of S2 on the linear chains is straightfor- 
ward, while it should be remarked at the junctions. The 
combination 1 — • tj in ^2 at the junctions is defined as 
follows: Three combinations are defined at the q = 6 junc- 
tions, and 5/2 combinations are defined at the q = 5 ver- 
tices. Figures[ljd) and[lje) show the combinations 1— ti-t2 
and 1— ti ■ (t2-l-t3)/2, which are respectively typical of the 
q = 6 vertices and of the q — 5 junctions. The combina- 
tion 1— ti • t2 is included in the summation of ^2 with the 
weight of 1, while the combination 1— ti ■ (t2-l-t3)/2 is 
included in ^2 with the weight of 1/2. This is the reason 
why the q = 5 junctions have 5/2 combinations. 
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The statistical mechanical model is defined 
tition function 

./ N 

l[dX^exp[-S{X)], 

where j' denotes that the multiple integration! 
performed under fixing the center of mass of t 
The integrations are simulated by the canon 
Carlo technique. The vertex position Xi is s' 
new position X^ with a three-dimensional ranc 
5Xi such that Xl = Xi + 5Xi . The random vei 
chosen in a small sphere, and the radius of th 
fixed at the beginning of the simulations so 
about a 50% acceptance rate; the radius of 
depends on b and on whether the surface is in 1 
phase or in the linear phase. 

Table 1. The surface size assumed in the simulal 
sizes {N,Ns,Nj) are assumed for each L' . 



L' 


{N,Ns,Nj) 


{N,Ns,Nj) 




8 


(2562,882,42) 


(5762,1982,92) 


(102 


6 


(3242,1442,92) 


(5762,2562,162) 


(901 


4 


(2562,1602,162) 


(5762,3602,362) 


(102 


3 


(3242,2522,362) 


(5762,4482,642) 


(90C 


2 


(2562,2562,642) 


(5762,5762,1442) 


(1024 



Table [T] shows the surface size {N, Ns, Nj 
in the simulations in this paper. Three differei 
assumed for each of the distance L' ranging f 
to L' = 8. 

The total number of MCS for the therm; 
2 X lO'' 3 X 10^ in the smooth spherical pi 
surfaces. To the contrary, the speed of converge 
low in the linear phase especially in the cases L' ='Z and 
L' = 3. In fact, the maximum number of the thermalization 
MCS in the linear phase close to the transition point is 
about 3 X 10^ for the surface of TV = 9002, L' = 3, and 
2 X 10^ for the surface of TV = 10242, L' = 2. In the cases of 
L' = 4 L' = 8, the thermalization MCS is about 6 x 10* ~ 
1 X 10^ in the linear phase close to the transition point 
on the largest size surfaces. Relatively smaller number of 
MCS is done for the thermalization in the smaller sized 
surfaces in each case of V . The total number of MCS for 
the production run is 2 x 10* ~ 3 x 10* on the smallest 
surfaces, and 3 x 10* 4 x 10* on the middle sized and 
the largest sized surfaces in each phase and in each L' . 



3 Results 

We firstly show in Figs. [2][a)-[2Ud) snapshots of surfaces 
in the smooth phase of size (a) N = 10242, L' = 8, (b) 
iV = 9002, L' = 6, (c) 7V = 9002, L' = 3, and (d) 7V = 10242, 
L' = 2. These snapshots are obtained at (a) 6=1760, (b) 



(a) (b) 




Fig. 2. (Color online) Snapshots of surfaces in the smooth 
phase obtained at (a) &=170, (b) & = 220, (c) & = 460, and (d) 
b = 480, where the surface size are (a) = 10242, L' = 8, (b) 
iV = 9002, L' = 6, (c) N = 9002, L' = 3, and (d) iV = 10242, 
L' = 2. Four snapshots were drawn in the same scale. 



6 = 220, (c) 6 = 460, and (d) 6 = 480. These snapshots were 
drawn in the same scale. From the surface sections, which 
are not depicted, we see that the surfaces in the smooth 
phase are empty inside the surface. Thus, we understand 
that the compartmentalized surface model in this paper 
has a well-defined smooth phase at sufficiently large 6 at 
least when L' <8. 

Figures[3ja)-[3fd) show snapshots in the linear phase at 
(a) 6=160, (b) 6 = 210, (c) 6 = 440, and (d) 6 = 460, which 
are close to the smooth phase in each case. The surface 
sizes (iV, L) in Figs. [3l^a)~[3l^d) are identical with those in 
Figs. [2l[a)-[2l[d) . The linearly extended surfaces are very 
different from each other in the length, so the scales of 
the figures are. We find that the linear phase is also well- 
defined. The planar phase observed in the model of [29], 
which corresponds to the model in this paper under the 
condition L' — l, is only seen on the surfaces of iV = 2562, 
L' = 2 and N = 3242, L' = 3. However, the planar phase 
does not appear on large sized surfaces even when L' = 2 
or L' — 3. We should note that there appear burs or pins 
on both ends of the surface in Fig. [3fc) although they are 
almost invisible. These burs or pins can only be seen on the 
surfaces of L' = 3 in the linear phase close to the transition 
point, however, the model is considered to be well-defined 
because of the finiteness of the physical quantities. 

The distribution of vertices and hence the surface shape 
is reflected in the mean square size X'^ defined by 

i i 



connectivity surfaces supported by skeletons 



(a) 



(b) 





Fig. 3. (Color online) Snapshots of surfaces in the linear phase 
obtained at (a) b = 160, N = 10242, L' = 8, (b) b = 210, iV = 
9002, L' = 6, (c) b = 440, = 9002, L' = 3, and (d) b = 460, 
A'^ = 10242, L' = 2, where the surface sizes are identical to 
those in Figs. [2ja)-{2jd), respectively. The scales of the figures 
are different from each other. 
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Fig. 4. The mean square size X vs. b obtained on the surfaces 
of size (a) L' = 8, (b) L' = 6, (c) L' = 4, (d) L' = 3, and (e) 
L' = 2. The solid symbols in (b) correspond to ill-defined data, 
which will be confirmed later in Fig. |8] 



where X is the center of mass of the surface. FiguresHJ^a)- 
mje) show X^ vs. b obtained on the surfaces of L' = 8, 
L' = 6, L' = 4, L' = 3, and L' = 2. We should note that 
the model has no collapsed phase because the model is ill- 
defined at sufficiently small b. The vertices are distributed 
extremely large volume in and behave just like a gas 
at We know that the model is well-defined in the 

whole region of b including 6^0 at L' = 1, and the well- 
definedness holds even at L' < 3. On the contrary, it is easy 
to see that the model turns to be ill-defined at L' > 4 at 



6-^0. The solid triangular symbols in Fig.lDJb) correspond 
to the ill-defined data obtained on the N = 5762 surface. 
This ill-dcfinedness is in fact seen in a large violation of 
Si/N — 3/2 and is confirmed later. 

The Hausdorff dimension H is defined by X^ ~ N^/^ 
in the limit of ^ oo. It is apparently expected that 
Hs ^ 2 in the smooth phase in all cases of L' from L' = 8 
to L' — 2. Then, we obtain in the linear phase by using 
X^ in that phase at the transition point. Table [2] shows 
Hs in the smooth phase and in the linear phase at the 
transition point. Hs in Table [2] is almost exactly identical 
to Hs = 2 the expected topological dimension of smooth 
surfaces, while Hl has values such that 1 < Hl < 2. We 
note that Hl should be Hl = 1 if the linear surface has 
constant density of vertices per unit length. The reason of 
the deviation of Hj^ from Hl = 1 in Table [2] is because the 
linear surface is not purely one-dimensional. 



Table 2. The Hausdorff dimensions Hs in the smooth phase 
and Hl in the linear phase. 



L' 


Hs 


Hl 


8 


2.00 ± 0.04 


1.57 ±0.02 


6 


2.00 ± 0.05 


1.54 ±0.11 


4 


2.01 ±0.04 


1.65 ±0.18 


3 


2.01 ±0.04 


1.50 ±0.03 


2 


2.00 ± 0.03 


1.75 ±0.09 




(c) \V=^ 
N=5762 




^(d)[r=3] 

N=5762T 

N=9002 to 



Fig. 5. (a) The one-dimensional bending energy S2/N's vs. b 
obtained on the surfaces of (a) L' = 8, (b) L' — 6, (c) H = 4, 
(d) L' = 3, and (e) L' = 2. 



Figures[Ha)-(S]^e) show the bending energy 5*2 /A^^ vs. b 
obtained on the surfaces of L' = 8 ~ L' = 2. The symbol Ng 
is given by N's =Ns + 2Nj — 6, which can also be written 
as Ng = 10^^ — 60m^ + 2. We see an expected jump in 
each S2/Ng at intermediate b, where X^ discontinuously 
changes. 
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Fig. 6. The two-dimensional bending energy Ss/Nb vs. b ob- 
tained on the surfaces of (a) L'=8, (b) L' = 6, (c) L' —4, (d) 
L' = 3, and (e) L' — 2. The solid symbols in (b) correspond to 
ill-defined data, which will be confirmed later in Fig. [S] 



The two-dimensional bending energy S3 is defined by 
53 = ^(l-n,.n,), (4) 

where is a unit normal vector of the triangle i, and 
y^^-^ denotes the summation over all nearest neighbor pairs 
of the triangles ij. S^/Nb vs. b is plotted in Figs. El^a)- 
inUe), where Nb is the total number of bonds including 
the bonds in the skeletons. Discontinuity in S^/Nb is very 
clear just as in the model of [29] , where L' = 1. We see that 
S^/Nb is very small in the smooth phase in the case of L' = 
2. This is because the surfaces are very smooth in that case 
just like the model in Surface fluctuations are almost 
completely suppressed in the smooth phase when L' = 2. 
This is expected from the compartment structure in the 
case of L'=:2 as shown in the snapshot of Fig.[T]Jc). In fact, 
there is no vertex inside the compartments on the surfaces 
of I/' = 2. On the contrary, the total number of vertices in 
a compartment is greater than one at least in the cases 
L'>3, and no bending energy is assumed on the triangles 
inside the compartments in those cases. This is the reason 
why Ss/Nb is relatively large in the smooth phase in the 
cases L' > 3. It should also be emphasized that the well- 
definedness of the surface is remarkable. Only constraint 
for the vertices and the triangles inside the compartments 
is the area energy 5*1 in our model, and no bending energy 
is assumed on those parts. The solid triangular symbols 
in Fig. mb) correspond to the ill-defined data obtained on 
the A^ = 5762 surface. 

We see that the transition point bc[kT] increases with 
decreasing L' . To see the dependence of physical quanti- 
ties on L' more clearly, we plot the bending energy S'2/iVg 
and the two-dimensional bending energy S^/Nb against 
L' in Figs. [7{a) and EJb), where the surface size is as- 
sumed as A^ = 5762 for all L' . We find from Fig. [Tja) that 
5*2 /-/Vg continuously changes against L' . On the contrary, 
the value of S^/Nb against L' in the smooth phase appears 
to be zero (non-zero) when L' < 2 {L' > 2). The reason of 
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Fig. 7. (a) The bending energy S2/N's vs. L', and (b) the two- 
dimensional bending energy Ss/Nb vs. L' . The surface size is 
iV = 5762. 



this behavior in S^/Nb is because the smoothness of the 
surface changes almost discontinuously depending on the 
compartment size L' as mentioned above. 
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Fig. 8. The Gaussian bond potential Si/N vs. b obtained on 
the surfaces of (a) L' = 8, (b) L' = 6, (c) L' = 4, (d) L' = 3, and 
(e) L' — 2. The solid symbols in (b) are obtained at 6 = 50 and 
6 = 60 and denote that the model is ill-defined at the rigion of 
6<60. 



Finally, the area energy Si/N is plotted against b in 
[5{a)-[HIe) in order to see the expected relation Si/N ~ 
3/2 is satisfied. If the surface generated by MC simula- 
tions is a well-defined two-dimensional object in R'^ then 
Si/N = 3/2 is expected from the scale invariant property 
of the partition function. We see in the figures that the 
relation is almost satisfied in the cases L' = 8 and L' = 6 
except the solid symbols, while the expected relation is 
slightly violated in the linear phase when L' < 3. We know 
the violation of the relation in the model at L' — > 1, and 
therefore the observed violation is consistent to the known 
result. The reason of a slight violation of the relation is, 
as described in , that the vertices move only in a one- 
dimensional region in R'^ in the linear phase, and hence 
the simulation for the three-dimensional integrations is 
not always well performed on such a linear surface. 

To the contrary, the solid triangular symbols on the 
N = 5762 and L' = 6 surface in Fig. [Hl^b) clearly represent 
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that the data are iU-defined because Si /N are far different 
from 3/2, where & = 50 and b — 60. The model is well- 
defined even at & = 40 on the N — 5762 and L' — 8 surface as 
we see from Fig.jSja), therefore, the lower bound b of well- 
deflnedness is not simply dependent only on L' . However, 
the model of L' > 4 is ill-defined when 6 in contrast to 
the cases L' < 3, where the model is well-defined even at 
6^0 as mentioned previously in this section. 



4 Summary and Conclusions 



In this paper, we introduced a compartmentalized struc- 
ture in a Nambu-Goto surface model on triangulated spher- 
ical surfaces. The mechanical strength of the surface is 
given by the one-dimensional bending energy defined only 
on the compartments with junctions. No bending energy 
is defined on the triangles inside the compartments, while 
the area energy term is defined all over the surface. The 
size of the compartment is characterized by the integer 
L' , which is the total number of bonds between the two 
neighboring junctions. We previously reported that a large 
variety of phases can be seen in a Nambu-Goto surface 
model, which is obtained from the model in this paper in 
the limit of L' ^ 1. In this sense, the model in this pa- 
per is an extension of the model in [55]. Therefore, it was 
non-trivial whether or not the compartmentalized struc- 
ture can be introduced into the model in ^29j , because the 
model might be ill-defined due to the compartments even 
at sufficiently large bending rigidity b. In fact, the model 
is obviously ill-defined in the limit of i'— >oo, where nei- 
ther the compartment nor the bending energy is defined 
on the surface. 

Five different sizes L' of compartments are assumed 
such that L' = 8, L' = 6, L' = 4, i' = 3, and L' = 2. For 
each L', three different sizes N of surfaces are assumed in 
the range A^ = 2562 ~ A^= 10242. By using the canonical 
MC simulation technique we study whether the compart- 
mentalized model is well-defined or not, and moreover we 
study the phase structure of the model at relatively large 
b region, where the smooth phase is expected if the model 
is well-defined. 

Our first observation is that the model is well-defined 
in the large b region if L' < 8 at least. In the cases of 
L' = 3 and L' — 2, the model is well-defined in the whole 
region of b, while it is not in the small b region when 
L' > 4. The second observation is that the model under- 
goes a first-order shape transformation transition between 
the smooth spherical phase and the linear phase at those 
relatively large b region. The planar phase, which is ob- 
served between the spherical phase and the linear phase 
in the model of L' = 1 in ^\ , disappears from the com- 
partmentalized model with L' >2. Moreover, no surface 
fluctuation is observed in the smooth surface in the case 
of L' = 2 just as in the model of |29j. The surface fluctu- 
ations of the smooth spherical surfaces emerge only when 
L'>3 and are completely suppressed when L' <2. 
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